Implementation of quantum search algorithm using classical Fourier optics 
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We report on an experiment on Grover's quantum search algorithm showing that classical waves 
can search a A'^-item database as efficiently as quantum mechanics can. The transverse beam profile 
of a short laser pulse is processed iteratively as the pulse bounces back and forth between two 
mirrors. We directly observe the sought item being found in ~ \/N iterations, in the form of a 
growing intensity peak on this profile. Although the lack of quantum entanglement limits the size 
of our database, our results show that entanglement is neither necessary for the algorithm itself, nor 
for its efficiency. 
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Quantum computers j^, ^ hold the promise of perform- 
ing tasks 1^, |4j that arc cither impossible or much less ef- 
ficient without the use of quantum mechanics. One such 
task is quantum searching introduced by Grover 

in a paper entitled " Quantum mechanics helps in search- 
ing for a needle in a haystacK^ . Consider using a phone 
book with N entries to find the name of a person whose 
phone number you have. Classically, this would require 
^ N consultations of the phone book. Grover's algo- 
rithm finds the desired entry with only ~ \/N consul- 
tations, using quantum mechanics. Here we show ex- 
perimentally that classical waves can find a needle in a 
haystack as efficiently as quantum mechanics can. Al- 
though some previous experiments ^, 0, ||, ||, |o| have 
demonstrated various aspects of quantum searching, all 
of them have been limited to four entries 0, |^ or a 
single query ||, ^, |T^. Our experiment closely follows 
Grover's algorithm, implementing for the first time an 
iterative search on a 32-item database, using classical 
waves. It provides a striking demonstration that the al- 
gorithm itself requires only wave properties |jll| but no 
entanglement Q. 

In Grover's (first) algorithm [Q each database item is 
associated with a quantum state. Initially the system is 
prepared in a superposition of all TV quantum states. The 
algorithm then amplifies the probability amplitude of the 
state being sought, in an iterative way. The item has 
been found once the probability amplitude of this "tar- 
get state" is near unity. Ideally this requires (7r/4)ViV 
iterations of the following two steps. In the first step a 
so-called "oracle" marks the item by inverting the phase 
of the associated quantum state [|| . In the second step 
the amplitudes of all states are inverted about the aver- 
age amplitude (lAA operation), converting phase infor- 
mation into amplitude information. 

The above protocol maps onto our classical-wave ex- 
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FIG. 1: Cavity implementing Grover's algorithm using opti- 
cal interference. We launch a short laser pulse with a gaussian 
transverse beam profile E{x), x representing the data regis- 
ter, into the cavity formed by mirrors Mi, 2. A line shaped 
depression in the oracle plate marks the item by imprinting 
a phase profile $o(a;). The sequence F^fFF^fF performs 
the inversion about average (lAA) as required by Grover's al- 
gorithm. Here F denotes a Fourier transform, performed by 
the lenses Li,2 (focal lengths fi = 400 mm, /2 = 600 mm). 
The lAA plate imprints a phase profile ^f{x') in the Fourier 
plane of the oracle. The enlargements show cuts of the phase 
plates perpendicular to the lines. As the pulse bounces back 
and forth, the transverse beam profile is processed iteratively 
and light is concentrated into the shaded mode. A high inten- 
sity peak, growing on the beam profile in the output plane, 
indicates the sought item. 



periment as follows (see Fig. |l|). A complex electric field 
amplitude E{x), viz. a transverse laser beam profile plays 
the role of the quantum probability ampliljudes. The con- 
tinuous coordinate x labels the items of the database, cor- 
responding to all possible quantum states. By spatial fil- 
tering we initialize the beam profile |i?(a;)p as a smooth, 
near-gaussian, distribution with a 1.33 mm diameter 
(FWHM; full width at half maximum). A single, ^ 300 ps 
laser pulse (wavelength 532 nm) enters a standing-wave 
cavity of 2.02 m optical path length through input mirror 
Ml (transmission 2%). The pulse travels back and forth 
between the cavity mirrors in 13.5 ns, each roundtrip rep- 
resenting one iteration of the search algorithm. Inside the 
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FIG. 2: Amplitude amplification as observed in trains of 
light pulses. Trace A shows the pulse train coupled out of a 
bare cavity, i.e. with the phase lines shifted out of the beam. 
The peak amplitudes decay exponentially with 25% roundtrip 
loss. For trace B the phase lines were moved into the beam 
and the pulses were recorded behind a narrow slit placed in 
the image of the oracle phase line. The energy in the pulses 
increases, even though the total energy decays. 



cavity an "oracle plate" marks the item by imprinting 
a phase profile on the beam, E{x) E{x) exp{i^o{x)), 
where $0(2;) — (j) in a narrow area around the "item po- 
sition" Xo and ^o{x) — elsewhere. Next, the lAA oper- 
ation is performed by the sequence F^fFF^fF, where 
F denotes a Fourier transform and denotes a phase 
plate like the oracle, but now imprinting a phase pro- 
file in the Fourier plane. The Fourier transforms 
replace the Walsh-Hadamard transforms ||T^ in the orig- 
inal proposal Q and are experimentally performed by 
spherical, achromatic doublet lenses Since F^ is a 
spatial inversion and $/(a;') = <&/(— x'), the lAA opera- 
tion reduces to F~^^^F. Thus the amplitude amplifying 
Grover iterator is ^lF~^^jF. Note that F^fF can be 
recognized as phase contrast imaging. 

We observe the progress of the search algorithm iter- 
ation by iteration, using the 2% transmission of mirror 
M2 after each cavity roundtrip. This light is imaged onto 
a 55 fj,m wide movable slit and the transmitted light is 
collected on a photodiode. The photodiode signal is am- 
plified and recorded by a digitizing oscilloscope. The 
light pulses are short compared to the roundtrip time, so 
that a train of output pulses is obtained, one pulse per 
iteration. In Fig. ^ we show two typical time traces. The 
trace in Fig. ^A has been recorded in an "empty cav- 
ity" , leaving the oracle and lAA plates inside the cavity, 
but moving the phase-shifting lines on the plates out of 
the beam. We observe an exponentially decaying peak 
amplitude, with a roundtrip loss of about 0.25, due to 
reflections. Next, we move the lAA phase line into the 
beam focus, put the oracle line in an arbitrary position 
in the beam, and place the detection slit in the image of 
the oracle line. We then observe a peak amplitude that 
grows during the first few iterations, even though the to- 
tal optical energy decreases. This is shown in Fig. I3.B 



and is a direct observation of amplitude amplification. 

We have measured the entire beam profile by recording 
traces like in Fig. ||.B for many different detection slit 
positions. We combined the peak values at the same time 
from different traces into a transverse beam profile. A 
sequence of such profiles for consecutive roundtrips shows 
how the algorithm proceeds. In Fig. 3.A-C we show 
three such sequences for increasing widths of the oracle 
line. Consecutive profiles within a sequence have been 
multiplied by a factor 0.75""'^, in order to compensate for 
optical losses. We clearly observe the solution growing 
as a high intensity peak in the transverse beam profile. 
The position of this peak is the position Xo of the sought 
item, i.e. the phase line in the oracle, imaged by the 
intracavity telescope. In the quantum case it would of 
course be impossible to watch the solution grow as the 
algorithm proceeds, because a measurement would cause 
the wave function to collapse. 

On the basis of Grover's algorithm we expec t the 
peak height to reach a maximum after {tt / 4) N / m 
roundtrips, where m is the number of marked items 
|]l5| , [1^ , and to oscillate through a sequence of maxima 
and minima with a period of (7r/2)-y/ N/m. In an ideal, 
loss-free system, these cycles of finding and "unfinding" 
would continue indefinitely. This period assumes that the 
phase shifts cj) have their ideal values. Since we use the 
plates in double pass inside the cavity, this ideal value 
is ±7r/2, whereas our measured value is (p — —1.1 ± 0.2 
rad. This increases the optimum number of iterations 
to [7r/(4sinl.l)],/7V/m. Although 4> may deviate from 
7r/2, a "phase matching" condition |17, 18) requires that 



the two phase shifts of the oracle and lAA plates are 
approximately equal. 

The ratio N/m can be interpreted as the size of the 
database for a single item search. Alternatively, the same 
N/m also describes a search for m adjacent items in a 
larger database of size N. The maximum database size 
is determined by optical diffraction, which limits the ef- 
fective number of positions x that can be resolved. For 
our cavity with a numerical aperture NA= 0.03, the 
limit on the resolution is given by Rayleigh's criterion 
as 0.61A/NA«10 ^m. For our 1.33 mm input beam, the 
maximum database size is then ~133. 

We can estimate N/m as the ratio of the input beam 
diameter to the oracle line width. The phase shifting 
lines have been produced as the shadows of thin metal 
wires (50, 100 and 200 fim diameter) while evaporating a 
thin layer of SiO onto a BK7 substrate. A phase-contrast 
image revealed line cross-sections that are well approx- 
imated by trapezoids, with flat inner regions of 42, 84 
and 126 /im, for the oracle plate and 136 /im for the 
lAA plate. The deviations are probably due to details 
of the evaporation procedure. Using the 1.33 mm diam- 
eter (FWHM) of the input beam, we get expected ratios 
N/m — 31.7, 15.8 and 10.6. We can compare this to 
the N/ m values as obtained from the position of the first 
maximum in the search, bearing in mind that the first 
image, having made 1/2 roundtrip, should be counted 
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FIG. 3: Iterative progress of the search algorithm as shown by measured and simulated beam profiles. Oracle lines of width 
(A) 42 /im, (B) 84 /im, and (C) 126 were used, corresponding to databases of 31.7, 15.8, and 10.6 items, respectively. Light 
coupled out of the cavity was recorded after each roundtrip through a scanning slit. The peak growing in the first iterations 
reveals the position of the sought item. The traces on the right (D, E, F) have been simulated using realistic experimental 
parameters, corresponding to those on the left. 



as 1/2 iteration. For the data shown in Fig. ||.A-C we 
estimate the maximum peak at 5, 3.5 and 3 iterations, 
leading to N /m = 32, 15.8 and 11.6 respectively, in good 
agreement wit h the expected numbers. The results thus 
confirm the N /m scaling behavior as expected from 
Grover's algorithm. 

The prime significance of the N/ni values is in the 
scaling of the searching period as y/N/rn. The absolute 
values of the expected N/ m depend on our chosen defini- 
tion for the input beam diameter (FWHM) and thus may 
seem somewhat arbitrary. Therefore we also compare our 
results to a simulation, shown in Fig. |3[D-F. We simu- 
late the phase plates by multiplying the beam profile by 
the experimentally determined trapezoidal phase profiles 
^o{x) and $/(a;). We describe the lenses by a Fourier 
transform. The results of the simulation agree well with 
the experiment, producing the maximum peak at the 
same number of iterations as the experiment. An im- 
portant difference between the experiment and the simu- 
lation is due to optical losses in the experiment. As men- 
tioned earlier, the experimental data have been scaled to 
compensate for the losses, which amplifies the noise in the 
last few iterations shown. Apart from this noise, we also 
see the development of side peaks. These are probably 
due to diffraction effects accumulating as the iterations 
progress, e.g. due to slight misalignments of our optical 
cavity. 

Keeping the resolution at ^10 ^m and extending the 




experiment to 2D, E{x,y), it should be feasible to per- 
form database searches of up to 10^ items experimentally, 
assuming a beam diameter of 1 cm. This is equivalent to 
about 20 qubits, so that we gain experimental access to 
problems that are as yet inaccessible for true quantum 
computers. These include quantum counting | [l6| , 
timation of the mean and median of a populationj 
the synthesis of arbitrary superposition states 
oretical studies have investigated fault tolerance 
and noise | |2^ in Grover's algorithm, predicting damp- 
ing of the cycles of finding and "unfinding" , like we also 
see in the experiment. The problem of "phase match- 
ing" ||l^, |8| can also be directly translated into optics as 
differential phase shifts provided by the oracle and lAA 
plate. These issues are as yet impossible to investigate ex- 
perimentally with present-day quantum computers. Our 
classical-wave experiment can bridge this gap. Note that 
it is complementary to a theoretical proposal by Farhi 
and Gutmann |Q to search a digital database in ana- 
logue time, rather than using discrete iterations. In our 
case, an analogue database is searched using discrete it- 
erations. 

Some classical- wave analogies of quantum information 
processing 25, p7|| , as well as a hybrid quantum- 
classical approach |2^ have been proposed previously. 
Some elements of Grover's algorithm have been demon- 
strated with classical waves The latter experiment 
demonstrated an oracle and lAA operation for a four- 
item database. Iterations were neither present nor nec- 
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essary, since for TV = 4 a single query reveals the sought 
item. A four-item database search has also been demon- 
strated using NMR techniques 0. Electronic wave 
packets in Rydberg atoms have been used to store and 
retrieve numbers B] and an equivalent experiment has 
been reported recently with classical light waves [^ . 
However, it has been pointed out that the Rydberg-atom 
experiment lacked the lAA operation p9| , which is a cru- 
cial ingredient of the quantum search algorithms. In our 
present experiment, Grover's second algorithm |^] can 
be recognized in the first transmitted pulse, which is es- 
sentially a phase-contrast image of the oracle. Since the 
contrast would be relatively low, the light pulse must con- 
tain sufficiently many photons to build up good readout 
statistics. By contrast, using Grover's first algorithm, 
the item could in principle be found with near certainty 
by sending a single photon through the oracle 0{\/N) 
times. 

It should be clear that our optical system is not a uni- 
versal quantum computer. Essentially we have mapped 
the 2"-dimensional Hilbert space of n qubits by the 
Hilbert space of a single photon in a superposition of 
2" transverse modes. It is well known [|ll|, ^ that 
this unary mapping comes at the cost of an exponential 
overhead in some physical resource. Previous classical 
analogies required an exponential number of components 
such as beam splitters |2^, |2^, The efficiency 

of a true quantum computer in implementing the trans- 
forms has been attributed to entanglement, i.e. to the 



tensor product structure of the Hilbert space. Despite 
the lack of entanglement in our present experiment, the 
Fourier transform is performed efhciently using only a 
single lens, independently of the size of the database. 
The lack of entanglement does however limit the size of 
the database, which scales linearly with the beam diam- 
eter D, or a for a 2D version. Thus the equivalent 
number of qubits scales only as oc log I?. Even if we 
set D equal to the size of the universe, ~ 10^^ m, this 
would yield only 206 equivalent qubits. This limitation 
exists for any database containing classical information. 
On the other hand, since Grover's algorithm provides 
only a ViV speedup, a quantum computer implementing 
Grover's algorithm becomes exponentially slow for an ex- 
ponentially large database. Thus our experiment shows 
that quantum entanglement is not needed to implement 
the algorithm or to improve the efficiency. Its only role 
in this case is to allow for a larger database size. 
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